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Abstract: In such cases it is necessary to construct exact quadratic formulas, as there are some errors 
in the approximate calculation of specific or non-specific integrals whose first-order derivatives 
become infinity in a given interval. In this work, the multiplicative method of loss of specificity under 
the integral of integral or non-specific integrals in which the first-order derivatives become infinity in a 
given interval is considered. 
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The problem of approximate integrals derived from a function with a finite number of products in the 
integration interval is well described in the mathematical literature. It is not possible to approximate the 
integral or non-specific integrals, one of the first derivatives of which is infinity, because it can give gross 
errors. Therefore, it is necessary to construct as accurate quadratic formulas as possible for such integrals. One 
such method is the multiplicative method, which allows you to get rid of the special under the integral. 


Here is the essence of the multiplicative method. 


Suppose, in the integration interval 


| o(x)dx (1) 


a 


Let the olx) function under the integral have some speciality. Suppose a function 


Q(x) = ox) f (x) 

be able to write in the form, where cx )is a positive function that contains all of the given properties, and 
f (x) is a (smooth) function with no special properties and 2n is an ordered product [2]. Using this, we can 
use Gauss's formula for the above integral (1) [3]: 


b 


J afx)f(x)ax= A(x)" (2) 


a 


while its residual hadi 
FOE Es 

R, (f) = Tay [Oe (x )dx, 

as&<b 


determined by the formula. Here, the function @ (x) is a polynomial of degree n, x is the root of ©, (x), 


n 


and A” is the coefficient of formula (2), which are: 
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œ, (x) 


1 f x 
(n) _ 
= oo} (x J ox) xf) 


The following properties are appropriate for these coefficients: 


P b 
l; Y A,” = | o(x)dx . 
k=1 f 


dx. (3) 


2. All A®” coefficients are positive. 


Thus, for the calculation of nonlinear integrals in a multiplicative way, all the properties of the integral are 
loaded into a weighted function, the corresponding orthogonal polynomial is restored, their nodes Pi are 


found, and Am are the coefficients[1]. 
Below we have a weighted function 


æ(x)=-—1nlxl, (-1<x< +l) 


We take the function and use the multiplicative method. 
In section [-1; 1], we construct polynomials @ (x), which are orthogonal to the weights 


ole - a i 


If we use the pair of the function cx) = —In|3} in [-1; 1], the polynomials @, (x) we are looking at are at 


(—1 <xs +1), and have the largest coefficients equal to one. 


level x“ , so that the rank indicator comes with n only in pairs. 


w (x) polynomials from polynomial theory 


n 


õp (x)= xð n (x)- 4n, (x) (4) 


n+l “n 


is found in the recurrent formula [1]. 
Of course, @(x)=1, &(x)=x become, &, (x)= x -4 from (4) 


Based on the property of orthogonality 

+1 +1 

| In|x|@, (x)dx =0 yoki [infall —A, Jax =0 
“| =i 

become, hence 


fin x(x? ~ A \dx =0 


l m ] 
Is derived. Let's split this integral down, Å, = 9° and then we find the @, (x) =x" - a That's the 


way we look at @(x). (4) based 
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From the condition of orthogonality 


[in — E +A, Je as =0 


from a pair of functions under the integral 


fina» -(5+4 a= 0. 
9 


Let's break it down and integrate it, 


> and, so, Õ (x)= x -< x. 


25 9 25 
That is, according to (4) 


9 l 

~ 4 2 

© (x)=x -| —+A, |x +—-/,. 
a)=at-[ 24 fe a 
From the condition of orthogonality 
+1 

[In|x|x°@, (x)dx =0, 

=] 

+1 


finx -x° Õõ, (x)dx=0, 


-1 


23-81 
T 
we define. In this way, the appearance of @, (x) polynomials 
H(x)=1, Gx)=x, Glx=x°-S, 


= o 4959-11657 , 25-2243 , 23-25-5689 
@,(x)= x FE, eee 
121.20347 11.20347 9-49 -121- 20347 
we define. Our next goal is to find the roots of the polynomials ®©, (x) found. To do this, using the 
Lobachevsky-Gref method, we find the roots of the polynomial @, (x). The coefficients of the constructed 


approximate formula are. 
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o, (x) 


: 1 
A! ) = ao x— x0” dx 


found using formula [1]. 


Table 1. 
Ea 
0 


p 2 | 0,3010300 _ 
a i 
a a 


a aa 
od 


aaa 1,4894550 
3 9 


o [11278121382716 | 0,1407330 


y” 


-0,748296| 0,135100 | 113903008 
, [0.224042 | 0,864900 1,9369658 


A” 
2 
I 
1 


0,221042 0,864900 1,9369658 


292 1351 — 
0,629296 0,135100 T.1303008 


-0,809432 0,064942 236125257 

-0,420490 0,387768 1,5885720 

0,000000 1,094580 0,0392476 

0,420490 0,387768 1,5885720 
432 4942 A 

0,80943 0,0649 3 8125257 


Based on the algorithm in the calculation, the calculation process of Table 1 was developed using high-level 
algorithmic languages (Delphi, Maple7), the roots of which are defined to 6 rooms. 


Conclusion: Based on the above algorithm, the process of calculating specific or non-specific integrals whose 
first-order derivatives become infinity in a given interval by multiplication is carried out without loss of 
specificity under the integral. 
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